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We consider a class of cosmological models in which the universe is filled with a (non-electric) 
charge density that repels itself by means of a force carried by a vector boson with a tiny mass. When 
the vector's mass depends upon other fields, the repulsive interaction gives rise to an electromagnetic 
barrier which prevents these fields from driving the mass to zero. This can modify the cosmology 
dramatically. We present a very simple realization of this idea in which the vector's mass arises 
. . . from a scalar field. The electromagnetic barrier prevents this field from rolling down its potential 

and thereby leads to accelerated expansion. 
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I. INTRODUCTION 



< 

^ - 

Recent observations of Type la supernovae with high redshift indicate that the universe is entering a phase of 
cosmological acceleration Identifying the causative agent is perhaps the most exciting task for fundamental 

theory at present. There are many candidates. It could be a cosmological constant, the need for which was suggested 
on the basis of other evidence even before the supernovae results |^-|| . Scalars will also work because one can construct 
a potential to support any homogeneous and isotropic geometry for which the Hubble constant does not increa se. pl 
Minimally coupled scalars becoming dominant at late times was also suggested before the supernovae results [R-ffoT 
Since then such models have been dubbed "quintessence" [jllj and have received extensive study |p^-|l6||. Nonminimal 
couplings have also been explored and recent inspiration has been derived from string theory p8| , p^ and extra 
■ dimensions |^,^ . It has even been suggested that quantum effects may be responsible ||2^ . 
', In the absence of compelling observational or theoretical support for any of the existing scenarios it is worth 
(~| considering what else might be driving the late time acceleration we seem to be seeing. We have long advocated that 



<^ I there might be interesting cosmological implications from long-range forces other than classical gravity [g3 24 1 . In the 
context of a universe which appears to be pushing itself apart, the repulsive force of vector boson exchange naturally 
^ ' comes to mind. That is, suppose the dark matter, or even ordinary matter, carries some charge — not electromagnetic 
f~| charge although we will henceforth use that language — that couples to a U{\) gauge boson which we shall refer to 

■ as "the photon." With an exactly massless photon the self-interaction would be infinite, but it can be made finite 
. ^ , by the simple device of endowing the photon with a tiny mass p5[ | . Since this mass can be much smaller than the 

■ Hubble constant it is still reasonable to speak of the resulting force as long-range. 

H Our purpose here is to survey the cosmology of such theories to see if they can lead to a late phase of acceleration. 

. . . 1 Therefore we shall not bother with the complications that would appear in a fully realistic or completely general 
model. However, it is important to note that there is no shortage of plausible candidates for the actual charge 
carrier, nor any lack of physical mechanisms for inducing the vector mass from a fundamentally invariant theory. One 
candidate for our charge density is dark matter itself because it is already believed to exist, although constraints from 
structure formation make such an identification problematic [ p6{ . A less restricted alternative, which we adopt here, 
is a minimally coupled charged scalar, similar to the field commonly postulated for "quintessence" models. The vector 
might acquire its mass, with a fundamentally U{\) invariant Lagrangian, through either spontaneous or dynamical 
symmetry breaking. The vector becomes massless when some order parameter vanishes. If we also enforce that 
increasing the order parameter increases the vacuum energy then potentially interesting cosmology can arise from the 
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tension between the repulsive interaction and the vacuum energy. It is neither necessary, nor even desirable at this 
stage, to commit to a specific model which illustrates this tension in a fully realistic fashion. 

In this paper we study the cosmological effects of such a long range repulsive force. We construct a simple model 
and analyze it, both analytically and numerically, within the framework of a homogeneous and isotropic spacetime. 
Our model consists of scalar QED with a nontrivial potential and an explicit C/(l)-breaking term which gives the 
photon a mass that depends upon the scalar's magnitude. Since the global U{1) is preserved, the associated charge is 
still conserved, and a homogeneous density of such charge experiences a repulsive self-interaction of the type described 
above. It turns out that the repulsive self-interaction cannot directly accelerate the universe, at least not for very long. 
However, it can do so indirectly through preventing the scalar from rolling down its potential and thereby driving the 
photon mass to zero. 

It should be noted that the simple model presented here is identical to "spintessence" ||2^,^ as long as only 
homogeneous and isotropic field configurations are considered. There are potentially important differences when one 
allows perturbations because our class of models contains a vector interaction which spintessence lacks. Note also that 
the model studied here is part of a much larger class. It is not necessary to imagine that either the charge density or 
the photon mass derive from complex scalars. We consider this mechanism only because it leads to a simple model 
in which accelerated expansion can be shown to occur under certain circumstances. 

This paper is organized as follows. In Section II we review how a slightly massive vector can be introduced without 
violating either homogeneity or isotropy. We also explain why the mass must depend upon other dynamical variables 
in order to obtain cosmologically interesting late time behavior. In Section III we study the simple model described 
above, first analytically and then numerically. This leads to the observation, not previously noticed for spintessence, 
that the most reasonable initial conditions result in a late time evolution in which the scalar bounces back and forth 
between its original potential and the electromagnetic barrier. We also discuss cosmological parameters. Section IV 
contains our summary and conclusions. 



II. VECTOR LONG RANGE FORCE 



A simple model for a repulsive cosmological force is a complex scalar field coupled to a U{1) gauge field. For 
example, one might choose a Lagrangian of the form 

g^''iD,^riD,^)~^g'^''gf'^F^pFp,-Vm) ■ (1) 

In a homogeneous fluid with a net charge density, the fluid will be self-repulsive. However, realizing this situation in 
an homogeneous, isotropic cosmology is problematic. For example, on a closed 3-manifold, the total charge of any 
infinite-range force field must be zero, so the charge density must vanish. (This is a consequence of the inability to 
define the "inside" and "outside" of an arbitrary Gaussian surface.) On an open manifold, it is possible to impose a 
nonzero charge density, but only at the expense of isotropy. One must choose a boundary condition at infinity which 
selects a direction for the lines of force. These obstacles can be evaded by simply making the vector massive |2^], for 
example with an explicit Proca term: 

C = Co + ^m^g'^-'A^A^V^ . (2) 

Note that the current density 

J^ = te[(j>iD^^)* -<j>*{D^^)] (3) 

is still conserved as a consequence of global U(l) invariance. 

The unique solution consistent with homogeneity and isotropy has to satisfy dicj) ~ and Ai — 0. Homogeneity also 
implies that diAo = 0, so that the field strength tensor i^^^ = dfj,Ai, — dvA^ vanishes. The only nontrivial equation 
of motion for the vector field is 




iecjfdacj) - ie(j)do(l)* + Aq {m^ + 2e'^(l)*(f)) = . (4) 
Therefore the unique solution for Aq is 

^ ie {(jydoct)* - 0*^00) , , 
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The question we wish to ask is: how does a spacetime dominated by such a charged scalar evolve? For simplicity, 
we choose a flat Friedmann-Robertson- Walker metric 



g^^i/dx^dx^ — dt^ — a?' [t) dx ■ dx . 
The evolution equation for the metric is the standard Friedmann equation 

Stt 



3TOpj 



Tf 



00 



(6) 



(7) 



The stress-energy tensor is defined by T^^y/^g = 25S/6g^'^ and its nonzero components in this geometry are given 
by the pressure and energy density, 



(8) 
(9) 



We see immediately that the repulsive interaction contributes a term which obeys p = p, rather than the p < — p/3 
needed for acceleration. This implies that the new term redshifts very rapidly, as a~^, and quickly becomes negligible. 
Physically this is because the Universe expands while the mass remains constant, so the force eventually becomes 
short range on cosmological scales. Therefore, the mere presence of a repulsive force does not generically lead to 
acceleration. It can do so only if the new interaction alters the scalar's evolution so as to make the potential 
more dominant than it would otherwise have been. 

We know that the self- interaction diverges in the massless limit. Imagine a situation in which the photon mass is 
not constant but rather depends upon some of the other fields. We might expect interesting cosmological effects if we 
begin with a nonzero charge density and then drive the vector mass towards zero. In this way the force can remain 
long range on cosmological scales even though the Universe keeps expanding. 

A simple model which allows this behavior can be obtained from the Lagrangian (0) by the replacement 2A^(/>^, 



g^^ {D^^r (D.^) + X^cj^^g^'A^A. - \g-^Pg''" F^pFp„ 



The homogeneous and isotropic solution (||) for the vector field becomes: 

ie{(t)do((>* - 0*9o0) 



2(A2 



The associated charge density is 



(10) 



(11) 



(12) 



Note that current conservation implies dQ{a^J'^) ~ 0. 

Homogeneity requires that the scalar depend upon time alone. With this simplification its equation of motion 
becomes 



3 - - 2ieAo 



dV 



(A2 + e2) Al + teAo + 3ie[^]Ao 



= 



Substituting (pj]) and making use of current conservation results in the form 



It is convenient to decompose the scalar into a magnitude and a phase, 

^(t)EE/(i)e»^W, 

so that Aq depends only upon the phase, 

e 



Ao 



A2 



(13) 



(14) 



(15) 



(16) 
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In these variables the real and imaginary parts of ( p^ become 

A2 



and 



/ + 3 



Of 



f 



A2 + e2 



2/ + 3 



-)/ 

a 



0, 



(17) 



(18) 



Up to coupling constants, the second equation is equivalent to current conservation and has the simple solution 

It is preferable to specify the initial charge density, which is conserved, rather than the value of the phase. The 
relation between the two variables is 



3 tO\ 



1 



A2 + e^ 



2eA2 J a^f 



Making this substitution for 9 the scalar equation of motion becomes 



/+3(-)/ 



K 



where the constant K is 



K : 



4e2A2 



(19) 



(20) 



(21) 



We see that the presence of the interaction gives rise to an "electromagnetic barrier" which prevents the scalar field 
from relaxing to the origin. 



The energy density and pressure are 



K 

<f2 

K 



Vif) , 
V{f) . 



(22) 
(23) 



Note that the interaction increases the energy density, as one would expect for a repulsive force. However, the 
interaction does not, by itself, act like "antigravity" in the popular sense of the phrase. The new term in the energy 
density and pressure looks like a fluid component with equation of state p = p, and redshifts like 1/a^. One might 
therefore expect that the interaction becomes negligible at late times. This is indeed true when the potential minimum 
occurs at some nonzero value of /. However, the situation is more interesting when the minimum is at / = because 
then the scalar is prevented from rolling down to its minimum by the electromagnetic interaction. In the next section, 
we discuss the dynamics of this case. 

III. DYNAMICS OF A MODEL WITH AN "ELECTROMAGNETIC BARRIER" 

In this section, we discuss the dynamics of a particularly interesting class of models with an "electromagnetic 
barrier" arising from a vector interaction of the type discussed in Section ^ Taking a Lagrangian of the form (p^), 
the equation of motion for the field / is given by Eq. (^0|). To complete the description of the dynamics, we must 
also include the Friedmann equation. 



Stt 
3m|[ 



K 



Vif) 



(24) 



We are particularly interested in potentials for which the minimum is at / = 0. For this reason, and for simplicity, 
we will assume a monomial form: 



V{f) = Vof' 



(25) 
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The full set of equations describing the dynamics of the system is 



-bVof 



h-l 



and 



Stt 
3mpj 



K 



(26) 



(27) 



It is convenient to scale the field f{t) so as to absorb the overall (87r/3m|[), and then redefine the couplings K and 
Vq to include the residual powers of this factor, 



If Stt 



\3mpj 



3m|j 

Stt 
3m|[ 



fit) 



K 



With these definitions, the equations of motion (Eq) and (ET^) read: 



and 



IbvF"-^ = 0, 



vF' 



3^2 



(28) 
(29) 
(30) 

(31) 
(32) 



Note that these equations are those of a field moving in an effective potential which depends upon the scale factor 
as well as the scalar magnitude, 



V,s [F, a) = vF' 



j6^2 ' 



(33) 



(This is similar to the time-dependent potential proposed for Variable Mass Particles or "VAMPs" |2^.) Due to the 
dependence upon a the minimum is changing as the Universe expands. 



^^0 



2k 
bv 



1/(6+2) 



-6/(b+2) 



(34) 



Based on this observation, we try to find a self-consistent solution to the equations of motion before proceeding with 
a numerical analysis. 

The strength of the electromagnetic barrier is controlled by the initial conditions, i.e. the initial charge density, 
//we assume that the electromagnetic barrier and the potential energy are each much larger than the kinetic terms, 
then they must approximately balance and the scalar field has to be approximately Fq. Under the same assumption 
we neglect the kinetic terms in the FRW equation (|3|), and obtain]^ 



k 



a^Fg 



(X a 



-6fc/(6+2) 



(35) 



The solution such that a{t = 0) = is 



^In this paper we use a subscript (e.g., -ffo, ao) to indicate the lowest-order solution, not, as is a frequently adopted convention 
in cosmology, the current value of a quantity. 
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aoit) = ao{to) 



to 



(fc+2)/3fc 



(36) 



This solution is indeed self-consistent at late times. While the potential energy scales like a~^''/'^^"'"^\ the kinetic energy 
drops much faster, like a~^. While the solution (Fo,ao) in Eq. ( ^ ) and Eq. (36) is a solution, and is self-consistent, 
it is not guaranteed that it is the solution, i.e. a late-time attractor for a range of initial conditions. To answer this 
question, we resort to numerical solution of the equations of motion. 

We have investigated a range of choices for the exponent b. We find that the solution (^^ for the time evolution 
of the scale factor is indeed approached at late times for a wide variety of initial conditions. Thus, we are able to 
model practically any kind of cosmological matter which is interacting via a vector force by an appropriate choice 
of b. For example, 6=2 corresponds to pressureless dust; 6 = 4 is radiation. However, of particular interest are scalar 
potentials that can provide an accelerated expansion. Acceleration occurs depending upon the sign of p + 3p. Using 
Eq. (p3) and the scalar's equation of motion (31) we find 



Stt 

3TOp[ 



(p + 3p) 



l^dt 



,dF_ 
dt 



2(1 ~b)V{F). 



(37) 



Therefore, acceleration is possible if 5 < ij^ Note that these potentials are very steep as F ^ 0, so the electromagnetic 
barrier is most effective for them. In what follows we will primarily concentrate on this class of potentials. Nevertheless, 
most of our observations are valid for an arbitrary b. 

It is important to note, however, that the solution ( p^ ) for Fq describes the behavior of the field only for very 
particular boundary conditions. If we assume the solution (jsj) for Fq the Friedmann equation becomes 



ao{to) 
aoito) 



V [Fo (to)] 
'b + 



[ao ito)f [Fo (to)]' 

b/(b+2) 

[ao\to)\ 



2k\ 
b^J 



where we have neglected the kinetic term. The initial time is to. On the other hand, using Eq. (pq), we have 



aoito) 



b + 2 



tn 



^aoito) J V 36 

These two equations imply that the initial value for the scale factor needs to b' 

2 i'b + 2\ ( 36 



(38) 



(39) 



ao (^o) 



vti 



b + 2 



2 /2^\ 6/(6+2)1 (^+2)/66 



bv J 



(40) 



Eq. (^) then implies that the initial value for the field is 

Fo ito) 



_i _2 / 2 \ /b + 2 
> 6 + 2/ V 36 



1/6 



(41) 



Note that the solution (|3J) with this boundary condition is an exact solution to the equation of motion (|3^), but only an 
approximate solution to the Friedmann equation. Numerical analysis indicates that if to is large, the solution (Fq, ao) 
is a stable solution of the full system of equations. For arbitrary initial conditions, however, the field does not smoothly 
follow the minimum Fo of the effective potential, but instead oscillates about the minimum. The oscillations do not 
damp with time. Despite the oscillation, however, the behavior (p6h of the scale factor is approximately maintained. 



^ Of course no potential with 6 < 1 could appear in the Lagrangian of a renormalizable quantum field theory. It might 
conceivably derive from quantum corrections to the effective potential. 

*We note that in a flat universe the value of a is arbitrary. We choose this particular scaling because it it convenient for 
understanding the role of boundary conditions and for numerical evaluation of the equations of motion. 
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FIG. 1. Field F vs. time for b = 1/2. The dashed line is the solution @ for Fq. 

Fig. 1 shows the field as a function of time for b — 1/2, fc = 1, and v — 0.01, oscillating about the minimum Fq. The 
boundary conditions are chosen such that the boundary conditions for the field are given by the solution (^l|), but 
the scale factor is set at 0.9 times its value in Eq. (pO[). To determine the oscillatory frequency around the minimum 
Fq, consider finearized fluctuations AF{t) and Aa around the self-consistent solution (Fo,ao), i.e. 



F{t) = Fo{t) + AF{t) 
a{t) = aoit) + Aa{t). 



(42) 
(43) 



The Hubble parameter can then be expanded as 



a 



- 1- 



Ha + Hi, 



ao 



1 dAa 
ao dt 



(44) 



where 



1 (dAa ^ 
Hi = — \ — HoAa 



d (Aa/ao) 



From 



ao \ dt 'J dt 

it follows that the fluctuations obey the following linear equations. 



(45) 



^AF + 3i7o|A^+(5 + 2) 



k 



2HqHi 



6- 



k Aa 



„6p4 



AF 



-6- 



k Aa 



(46) 



where iJo(i)> oo(i)) -fo(0 the known unperturbed quantities. The Fq term is included in the expression because 
it is the perturbation which prevents the self-consistent solution (Fo,ao) firom solving the Friedmann equation. This 
term proves in practice to be negligible. 

The friction term can be eliminated through the rescaling 



AF{t) = g(i)«o 

The inhomogeneous equation (|4q) for the scalar field AF{t) reduces to 



(47) 
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dVrf^^+-^ff = -6af — , (48) 



where the square of the frequency is 

2 3 

" =4 



„dHo 2 



+ (& + 2)-6t;4- (49) 



The frequency therefore depends entirely on unperturbed quantities. It is straightforward to show that the first two 
terms are subdominant for b < 2, with 

77^ oc i/o oc ^'/(''+'\ (50) 

so that in the late time limit, 

7 /07 \ -4/(6+2) 

.2 ^ iu , o\ ^ -U(U,0\(^'^\ 6{2-b)/(2+b) 



We therefore expect the solution to be of the form of a driven oscillator with frequency given by Eq. (|l]) . Numerical 
solution agrees with this analytical result as shown in Fig. 2. 

The time dependence of the driving term on the right-hand side of Eq. ( ^ ) can be evaluated as follows. We take 
the perturbed Friedmann equation ( |46| ) and note that the term 

(X a-' (52) 

becomes negligible at late times. We then have 

c?(Aa/ao) 6fc / Aa 

II = z_nc 

Since we have the unperturbed solution 



we can write the Friedmann equation as 



We then have a solution for Aa, 



ao 

The time dependence of the oscillator equation (MS) can then be evaluated. 



\ o,a-''/^''+'Ut-\ (56) 



,2„_ „3/2 



9 + ^9^~ao 3 



6k / Aa 



0^0 \ ao 

3(6-56)/2(6+2) 



(57) 



The right-hand side of Eq. (^) therefore grows with time and drives the oscillations of the scalar field F. While we 
do not have an analytical solution to Eq. (|48|), we observe numerically that at late times AF/Fq ^ const., consistent 
with the rapid growth in amplitude of the oscillations of 5 = a^^^AF. Fig. 2 shows the numerical field solution as a 
function of the integrated frequency, in agreement with the analytical result. In fact, the driving of the oscillations 
and the frequency derived through perturbative analysis are quite robust even in the limit where the perturbative 
expressions above are no longer accurate. 
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FIG. 2. Field F{t) normalized to the analytic solution Fo{t) ^ vs. / Lj{t)dt, for b = 1/2. 

The most natural initial state is one in which the photon has a large and approximately time independent mass. 
In this model, this corresponds to the scalar having a large potential energy compared to its kinetic energy and 
electromagnetic barrier. This particular initial condition leads to an oscillatory late time behavior. The details of the 
oscillations (i.e. the average value of F) depend on the initial conditions, but the oscillations are generic. The field 
begins to oscillate even if it starts at the instantaneous minimum of the effective potential, unless the initial condition 
for the scale factor is given exactly by Eq. (^0|). The reason is that the time dependence of the effective potential 
serves as a driving term for the oscillations. Fig. 3 shows the three components of the stress-energy (potential, kinetic, 
and barrier) as the field evolves. Although the potential dominates throughout the evolution, all three terms scale 
identically with a. The precise relationship between the three components is very sensitive to initial conditions. Note 
in particular the fact that the kinetic energy scales identically with the other terms in the stress-energy, unlike in the 
case of the self-consistent solution Fq. 
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FIG. 3. The three components of the energy density vs. scale factor for b — 1/2. The potential (upper line, blue) and 
"barrier" (middle line, red) terms dominate the energy density, while the kinetic term (lower line, black) is subdominant. All 
three scale identically with a. 

Of most interest are, of course, the cosmological parameters. Do the small oscillations of the field affect the 
background? The deceleration parameter q follows the behavior of the field (Fig. 4). At times when the field is at 
its maximum, q reaches its minimum at —0.7. Note in particular that despite the fact that the kinetic energy of the 
field vanishes when q is at its minimum, the equation of state never reaches p = —p due to the contribution of the 
barrier term to the stress-energy, and q > —1 at all times. The system, however is on average in a state of accelerated 
expansion. 




ln(t) 

FIG. 4. The deceleration parameter q vs. time. The time variable is in the internal units of the numerical solution, 
t = [20, 100] 
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The behavior of the Hubble constant for the same set of initial conditions is shown in Figure 5. It, too, exhibits 
oscillations. Although the specifics of the field evolution are model dependent, they arise from the field oscillating about 
the minimum of the effective potential created by the potential and barrier terms. So, oscillations can be expected in 
any model of this type. Therefore the anomalous jitters in the Hubble constant and deceleration parameter can serve 
as a distinguishing feature of such models. A possible observational signature of such behavior would be the presence 
of scatter on a Hubble diagram which cannot be attributed to statistical error. 




3 3.S 4 4.5 



ln(t) 

FIG. 5. The Hubble parameter vs. time. The (red) dotted line is the result for the self-consistent solution. 

The Hubble parameter itself is not directly observable, however. Cosmological expansion is measured by comparing 
the redshifts and distances of standard candles such as Type la supernovae. The luminosity distance measured by 
the observed brightness of a known standard candle, is 

d^=a{to){l + z)£ (58) 

where to is the current time, and z is the redshift of the source. This is most commonly expressed as the distance 
modulus, or dimming (in magnitudes) of the observed object: 

(m-M) = 51og(dL/10pc) = 51og(iJodL) +41.6. (59) 

Here M is the absolute magnitude of the standard candle, m is the observed magnitude, and the luminosity distance is 
expressed in units of the Hubble constant today, H{to)di,. We have taken H{tQ) = 70 km/sec/Mpc. Since we consider 
a toy model in which only the oscillating scalar contributes to the energy density of the universe, the choice of what 
to label the current time to is arbitrary. (A more realistic model would consist of both matter and scalar field, and 
the current time would be fixed by the observed ratios of their densities, ftu ~ 0.3, fl^ ~ 0.7.) We take an optimistic 
choice for to, such that the early, large oscillations in the Hubble parameter are at a redshift of approximately 1 (which 
in the units used in the numerical simulation, is to = 68.) Fig. 6 shows a plot of distance modulus vs. redshift for the 
Hubble parameter plotted in Fig. 5. Despite the strong oscillations in the deceleration parameter, the oscillations on 
the Hubble diagram are very small, in fact too small to be evident on the plot. 
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FIG. 6. Distance modulus vs. redshift. Note that the oscillations in the Hubble parameter are not visible in this graph. 

The magnitude of the oscillations on the Hubble diagram can be shown by plotting the residual A(m — M) between 
the exact solution and the self-consistent solution, shown in Fig. 7. Fig. 8 shows a plot of the deceleration parameter 
g as a function of redshift for the same choice of to- 



< 




FIG. 7. Residual A(m — M) vs. redshift, relative to the self-consistent solution. The oscillations are visible at the level of 
A(m - M) ~ 0.04. 
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FIG. 8. Deceleration parameter q vs. redshift, for reference to Figs. 6 and 7. 

We see that the magnitude of the oscillations is A(m — M) ~ 0.04, smaller than the accuracy of current supernova 
la observations, which have errors on the order of A(m — M) ^ 0.15 — 0.20 js^. Other choices of initial conditions 
result in a similar order of magnitude for A(rn — M), to within a factor of a few. Detection of such a signature would 
present a formidable observational challenge, but it is conceivable that future measurements such as those by SNAP 
pl| could approach this level of accuracy. We note with interest that there is in fact some evidence of small scatter 
beyond observational uncertainties in existing data, at the level of .12 magnitude [Eol. 



IV. SUMMARY AND CONCLUSION 



In this paper we studied how a presence of a long range repulsive force affects cosmology. We have explicitly 
constructed a model for a cosmology dominated by a charged scalar field with a long-range repulsive interaction. The 
simplest such model, one with an explicitly broken U{1) gauge interaction (|^) is sufficient to show that a universal 
repulsive interaction does not generically introduce a negative pressure, but rather a term with equation of state 
p = p. Thus, a universal repulsive force by itself does not act like a cosmological constant. It can have an interesting 
effect on cosmology, including support of inflation, but in a more subtle way. 

In order for the repulsive force to be of any significance at late times, its range has to be allowed to increase, i.e. 
the mass of the force carrier needs to be dynamical. Then because of charge conservation, an electromagnetic barrier 
arises |2^ (in addition to a kinetic barrier as in spintessence [ p7|j28| ) . The system behaves as a driven oscillator around 
the effective potential arising from the sum of its original potential and the electromagnetic barrier. Even though 
the resulting electromagnetic barrier contributes a. w = +1 term to the stress tensor it can lower the deceleration 
parameter q by pushing the field back up its potential. For concave potentials, the deceleration parameter can even 
be negative. 

Here we have chosen the simplest case of a dynamically changing mass, that is, one without any additional fields 
required. Because of this choice, the electromagnetic barrier has the same form as the kinetic energy barrier and 
they can be combined. The resultant equations of motion thus look like spintessence. An important distinction arises 
when spatial perturbations are considered. In our model this would generate a nonzero electric and magnetic fields 
which would tend to resist clumping of the charges. Therefore, it might distinguish our model from these models with 
respect to their tendency to decay into Q-balls Q . 

We study the model both analytically and numerically. Our analytical formulas are confirmed by the numerical 
simulations. We find that there exists a self-consistent solution such that kinetic energy is negligible compared to 
the potential energy of the scalar field and its electromagnetic barrier (identical, in this particular model, to the 
spintessence solution). However, for generic initial conditions, the solution displays oscillatory behavior. We find that 
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the small oscillations change scaling of the kinetic energy, so that it scales in the same way as the potential and barrier 
terms in the stress-energy. Depending on initial conditions, the kinetic energy can be of the same order of magnitude 
as the other terms in the stress energy, unlike the self-consistent solution. 

The presence of oscillations can result in interesting observable consequences. In particular, the Hubble parameter 
and the deceleration parameter will also in general oscillate, creating an anomalous "spread" in a Hubble diagram 
not explainable by statistical error. Depending upon how big the effect is, it may or may not be observable and it 
may possibly constitute a signature of a long range repulsive force. 

In conclusion, we have presented the simplest of a class of models with a long range vector force. We emphasize 
that the key feature that makes our model work is the electromagnetic potential becoming increasingly long range. 
This is achieved by driving the mass of the vector dynamically to zero. However, apart from the simple case presented 
here, the same effect can be achieved through a wide range of possibilities. For example, we can have another scalar 
field generate the photon mass. We can have it generated by dynamical symmetry breaking of a Fermi field. We can 
even have the charge density provided by a Fermi field condensate. Nevertheless, many of the features present in this 
specific model, including the anomalous jitters, can be expected to be generic for models of this type. 
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